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Abstract: We present the "endpoint" formalism for the calculation of electromagnetic radiation and illustrate its applica- 
tions in astroparticle physics. We use the formalism to explain why the coherent radiation from the Askaryan effect is not 
in general Cherenkov radiation, as the emission directly results from the time- variation of the net charge in the particle 
shower. Secondly, we illustrate how the formalism has been applied in the air shower radio emission code REAS3 to 
unify the microscopic and macroscopic views of radio emission from extensive air showers. Indeed, the formalism is 
completely universal and particularly well-suited for implementation in Monte Carlo codes in the time- and frequency- 
domains. It easily reproduces well-known "classical mechanisms" such as synchrotron radiation, Vavilov-Cherenkov 
radiation and transition radiation in the adequate limits, but has the advantage that it continues to work in realistic, com- 
plex situations, where the "classical mechanisms" tend to no longer apply and adhering to them can result in misleading 
interpretations. 
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1 Introduction 

Processes leading to the emission of electromagnetic radia- 
tion can be very complex. To get a grasp of these processes, 
physicists often revert to "classical named radiation mech- 
anisms" such as synchrotron radiation, transition radiation 
or Vavilov-Cherenkov radiation and try to apply these con- 
cepts to the problem at hand. However, it turns out that 
this strategy can often yield misleading results. Usually, 
the classical mechanisms only apply to idealized problems 
involving, e.g., infinite particle tracks, and their application 
to complex realistic situations can cloud the view for the 
complete picture. 

Here, we propose a formalism suitable for the calculation 
of electromagnetic radiation from any kind of particle ac- 
celeration, which lends itself very well to the implementa- 
tion in Monte Carlo codes, and necessitates no simplifying 
approximations to the underlying processes. In this "end- 
point formalism", presented in detail in Ref. [ 1 1, the trajec- 
tory of individual particles is described as a series of points 
at which particles are accelerated instantaneously, leading 
to discrete radiation contributions which can then be easily 
superposed. 

The endpoint formalism can be applied in both the fre- 
quency and time domains, and correctly reproduces the 
aforementioned classical named processes. Furthermore, 



we demonstrate how the application of the endpoint for- 
malism has provided essential insights for radio emission 
processes in astroparticle physics. 

2 The endpoint formalism 

The endpoint formalism relies on the fact that any trajec- 
tory of a charged particle undergoing acceleration can be 
decomposed into a series of discrete "acceleration events" 
joined by straight tracks. The error introduced by this pro- 
cedure can be made arbitrarily small by applying the de- 
composition on scales as small as necessary. 

Furthermore, each of the acceleration events during which 
the velocity vector of a charged particle changes, can again 
be decomposed into two "endpoints", one corresponding 
to the instantaneous acceleration of the particle from its ve- 
locity to rest and the second corresponding to the instanta- 
neous acceleration of the particle from rest to its new ve- 
locity. 

The electromagnetic radiation associated with these end- 
points can be calculated with a simple formula, which 
lends itself well to application in Monte Carlo simulations. 
(For derivations of the stated equations we kindly refer the 
reader to (TJ.) 
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2.1 Frequency domain 

The electric field seen by an observer at position x and ob- 
serving frequency v from an individual endpoint is given 
by: 



E±(x,is) = ± 



q e 



ikR(t' ) 



c R(t' ) i - n (3* ■ f 
■ fx [fx 0*] (1) 

where k — 2-Kvnjc, q denotes the particle charge, n is 
the medium refractive index, f gives the line of sight from 
observer to particle, R = \ f\ and t' corresponds to the (re- 
tarded) time at which the endpoint radiates. f3* denotes the 
non-zero velocity associated with the endpoint. '±' is pos- 
itive when the acceleration is from rest to j3* and negative 
when the acceleration is from (3* to rest. 

2.2 Time domain 

Likewise, the electric field an observer sees from an end- 
point in the time domain is given by: 



E±(x,t) = ± 



1 q I f x [f x j3* 
Ai7 



(l-n/3* -f)R 



(2) 



Here, the result has to be interpreted as the electric field 
time-averaged over a time scale At, the adequate choice of 
which is dictated by the time resolution of interest. 

2.3 Example application 

In Fig. Q~|we illustrate the electromagnetic radiation asso- 
ciated with the acceleration of individual charged particles 
as calculated using the endpoint formalism in the frequency 
domain. In the figure, the emission from processes such as 
the instantaneous acceleration of a particle from rest ('ac- 
celeration'), the deflection of a relativistic electron by 20° 
('deflection'), the deceleration of a relativistic electron to 
rest ('slow-down'), and the reversal of direction of a mildly 
relativistic electron ('reversal'), is shown. Typical effects 
such as the beaming of radiation for relativistically moving 
charges and the shock wave behaviour near the Cherenkov 
angle occurring in media with a refractive index different 
from unity turn up naturally in this calculation. 

3 Application to classical mechanisms 

To demonstrate that the endpoint formalism can reproduce 
"classical named radiation mechanisms", we briefly illus- 
trate its application to the processes of synchrotron radia- 
tion and Vavilov-Cherenkov radiation. 

3.1 Synchrotron radiation 

While the term "synchrotron radiation" usually implies 
relativistic particles revolving continuously in a magnetic 



1e-15 



1e-16 



vacuum (n=1) 



£ 1e-17 



1e-18 



E 
> 



1e-1S 



1e-13 



1&-14 



1e-15 



Acceleration to y= 100 
Deflection through 20° (y = 10) 
Slow-down (7=11 to 1 0) 
Reversal (y = 2) 



30 60 90 120 150 

Angle 6 from initial velocity vector (degrees) 



180 



£ 1e-16 



1e-17 : 



1e-18 



1e-19 



dielectric (n=2) 


Acceleration to y= 100 


Deflection through 20° (y= 10) 




Slow-down (7 = 11 to 10) 




Velocity Reversal (y = 2) 


* „c * 


A *cA 


o 
t I 

n : 
ii :■ 
[i :: J 
' 1 I'.l 





1 1 ; 

' \1>t\ 





30 60 90 120 150 

Angle 8 from initial velocity vector (degrees) 



180 



Figure 1 : Electric field magnitude resulting from the accel- 
eration of a relativistic electron in four simple cases (see 
text) in (top) vacuum and (bottom) a dielectric with refrac- 
tive index n = 2, as calculated using Eq. ([TJ. 



field, here we look at a particle undergoing precisely one 
revolution. To describe this situation using the endpoint 
formalism, the helical particle track is decomposed into a 
series of straight tracks joined by "kinks" in which the par- 
ticle is accelerated instantaneously. As explained earlier, 
each kink in turn is represented by two endpoints with con- 
tributions and E + , describing the instantaneous decel- 
eration to rest followed by the instantaneous acceleration 
to the new velocity (in this case with the same magnitude, 
but different direction). 

Using either of equations ([T} and (fJJ, the emission radi- 
ated during the particle revolution can then be calculated 
as the superposition of the individual endpoint contribu- 
tions. Care has to be taken to describe the particle motion 
on a fine-enough spatial scale, as dictated (respectively) by 
the maximum frequency/minimum time interval of inter- 
est. The results for both the frequency- and time-domain 
are shown in Fig. [2] Note that calculating the result for one 
domain by an FFT of the result of the other domain yields 
the same result, but is usually much more involved than 
performing a direct endpoint calculation. 

Comparing these results with the established theory of syn- 
chrotron radiation confirms that the endpoint formalism 
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Figure 2: Time-trace of a synchrotron pulse produced by a single gyration of an electron with ft = 0.999 and gyration 
radius r = 100 m (top), and power frequency spectrum of the emitted radiation (bottom). In each case, the direct cal- 
culations (time- and frequency-domain respectively) are shown in comparison with the results generated by fast-Fourier- 
transforming data from the other domain. 



can reproduce the well-known theory. (For further details 
please refer to Ref. (TJ.) 

3.2 Vavilov-Cherenkov radiation 

In the scientific literature, the term "(Vavilov)-Cherenkov 
radiation" is often used somewhat imprecisely, but typi- 
cally refers to the original calculation of Frank and Tamm 
in the case of charged particles propagating through a di- 
electric medium along infinite tracks [2 1. There is no accel- 
eration associated with this particle movement, and there- 
fore this "true" Vavilov-Cherenkov emission cannot be de- 
scribed with the endpoint formalism discussed here. 

However, in realistic situations particle tracks are not infi- 
nite. When trying to calculate the emission associated with 
particles propagating through a medium on finite tracks, the 
emission is usually calculated using "track segments", and, 
as first calculated by Tamm [3 1, the far- field result near the 
Cherenkov angle turns out to be proportionate to the total 
length of track traversed. 

Using the endpoint formalism, we regard the radiation as 
contributions from the endpoints of the track segments. 
Calculating the radiation based on equation (Q}, we can de- 
rive (cf. Ref. [ T| for details): 



£track(£, v) 



q/3 sm e i(fefi+2™ti) 
c R 

p27rzi/(l — nj3 cos 6)St 

1 — nf3 COS 6 



E (3) 



where 9 denotes the angle between the particle velocity 
vector and the line of sight from observer to endpoint, t\ 
denotes the (retarded) time of the start of the track, 5t is the 
track length divided by c, and an approximation was made 
for observer distances large with respect to the track length. 



This result corresponds (within a factor of 2 arising from 
a different definition of the Fourier transform) to the 
"Vavilov-Cherenkov radiation formula" of Eq. 12 in Zas, 
Halzen, and Stanev [4|, with \i r = 1 and q = — e. These 
authors also show that their result — and, by extension, 
ours — produces the same emitted power as Tamm's cal- 
culation. Indeed, it has previously been shown that the ra- 
diation from a finite particle track consists of both endpoint 
contributions and a "true" Vavilov-Cherenkov contribution 
from the track itself, the latter of which is negligible at most 
angles in the far- field 1 5 1 . 

In other words, although the ZHS code is commonly un- 
derstood as calculating "the Vavilov-Cherenkov radiation" 
from a particle cascade in a dense medium, it actually cal- 
culates the radiation from the endpoints of finite particle 
tracks. As the particle cascade exhibits a variation of the 
negative charge excess during its evolution, the endpoint 
contributions of the individual tracks do not cancel and 
a net contribution arises. This radiation due to the time- 
variation of the net charge is a central, yet not generally 
appreciated aspect of the Askaryan radiation |6| process of 
cascades in dense media. Thus, we expect that the most re- 
cent results on "Cerenkov radio pulses" from straight track 
segments [7 |, wherein the authors use a method based on 
calculating the vector potential, is generally applicable to 
all radiation processes. 

4 Applications in astroparticle physics 

We have developed the endpoint formalism in the context 
of astroparticle physics and hereby illustrate how it has 
helped us to achieve a better understanding of radio emis- 
sion processes relevant for cosmic ray and neutrino detec- 
tion. 
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Figure 3: A radio pulse simulated with REAS3 (solid 
red line) compared with the radio pulse simulated with 
REAS2.59 (blue dashed line). REAS3 is based on the end- 
point formalism. 

4.1 Radio emission from air showers 

Radio emission from extensive air showers has been under 
intense study over the recent decade. To interpret the data 
gathered by experiments, a good understanding of the un- 
derlying emission mechanism is vital. A significant num- 
ber of models have been developed in recent years, but dif- 
ferent models gave different results. In particular, there was 
a discrepancy between microscopic Monte Carlo based ap- 
proaches such as REAS2 [ 8 1 and the macroscopic MGMR 
model |9|. It became clear that in the microscopic time- 
domain treatments, a radiation contribution associated with 
the time-variation of the number of charged particles had 
been neglected. Using the endpoint formalism in the time- 
domain in RE AS 3 [10], it could be understood in detail 
how the superposition of the individual particle emissions 
can be carried out consistently. Consequently, the differ- 
ent models could be reconciled IfTTI — a breakthrough in 
the modelling of radio emission from extensive air show- 
ers. The importance of the correction introduced by the 
endpoint formalism is illustrated in Fig. [3] 

4.2 Radio emission in the lunar regolith 

Numerous projects aim at detecting cosmic rays or neutri- 
nos via radio emission radiated from showers initiated in 
the lunar regolith. For a long time it has been the general 
understanding that the radio emission would vanish if the 
particle shower occurred close to the surface of the moon 
due to the so-called "formation zone effect" (first investi- 
gated in this context by Ref. El). The argument is that the 
refractive index close to the surface effectively approaches 
unity and that therefore Vavilov-Cherenkov emission has to 
vanish. 

As has been argued above, however, the emission pro- 
duced by showers in dense media (the Askaryan radiation) 
is dominated by radiation from the endpoints of individual 



particle tracks. The fact that the shower grows and declines 
again means that there is a net radiation contribution due 
to the variation of the net charge, which leads to radiation 
irrespective of the refractive index of the medium. A sup- 
pression of the emission due to the presence of the surface 
still occurs due to refractive effects at the surface — how- 
ever, this effect applies equally to both shallow and deep 
cascades, and is taken into account in experimental simula- 
tions. 

5 Conclusions 

We have developed an endpoint formalism for the calcu- 
lation of electromagnetic radiation in arbitrarily complex 
radiation processes. The formalism is universally applica- 
ble and can be implemented easily in calculations both in 
the time and frequency domains. 

Classical named radiation processes such as synchrotron 
radiation or Vavilov-Cherenkov radiation can be success- 
fully reproduced using the endpoint formalism. However, 
unlike such descriptions, the endpoint formalism does not 
make any simplifying assumptions which often break down 
in realistic situations. 

We have successfully used the endpoint formalism to ad- 
dress problems in astroparticle physics. In particular, we 
were able to calculate the radio emission from extensive air 
showers using the endpoint formalism, and thereby recon- 
cile discrepancies existing between earlier modelling ap- 
proaches. Also, the endpoint formalism illustrates that 
Askaryan radio emission in the lunar regolith does not van- 
ish even if the particle cascade develops close to the sur- 
face, as the net growth and decline of the shower will pro- 
duce radio emission even for a refractive index approaching 
unity. 

References 

[I] C. W. James, H. Falcke, T. Huege, M. Ludwig, Phys. Rev. E 
(submitted), larXiv: 1007.41461 

[2] I.M. Frank and I.E. Tamm, Dokl. Akad. Nauk SSSR 14, 109 
(1937). 

[3] I.E. Tamm, J.Phys.Moscow 1, 439 (1939). 

[4] E. Zas, F. Halzen, and T. Stanev, Phys. Rev. D 45, 362 (1992). 

[5] G.N. Afanasiev, V.G. Kartavenko, and Y.P. Stepanovsky, 

J.Phys.D 32, 2029 (1999). 
[6] G.A. Askaryan, Sov. Phys. JETP, 14, 441 (1962); 48, 988 

(1965). 

[7] J. Alvarez-Mufiiz, A. Romero- Wolf, E. Zas, Phys.Rev. D 81, 
123009 (2010). 

[8] T. Huege, R. Ulrich, R. Engel, Astrop. Phys., 2007, 27: 392- 
405 

[9] O. Scholten, K. Werner, F. Rusydi, Astrop. Phys., 2008, 29: 
94-103 

[10] M. Ludwig, T.Huege, Astrop. Phys., 2010, 34: 438-446 

[II] T. Huege, M. Ludwig, O. Scholten and K.D. deVries, 
ARENA Nantes, France (2010), Nucl. Instr. Meth. A in press, 
doi:10.1016/j.nima.2010.1 1.041 

[12] PW. Gorham, K.M. Kiewer, C.J. Naudet, D.P Saltzberg, 
and D.R. Williams, RADHEP 2000, |arXiv:astro-p h/0 1 02435 
(2001). 



